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1. INTRODUCTION 

The development of induction motor drives has a long history for decades. Many control algorithms 
have been developed to meet the requirements of the industry [1]. Besides, high power density has become a 
research trend in recent years in power electronics field. Accordingly, the new converters are designed not 
only to meet the power requirement but also archive a low volume and lightweight in many applications [2]. 
The development of semiconductor technology using GaN material allows them to operate at switching 
frequency from 50 to 200 kHz instead of only 2 to 20 kHz on IGBT transistors [3]. In addition, GaN-based 
devices also have advantages such as reducing switching losses and high temperatures operation. Many 
studies have been published to prove its effectiveness in practice with power electronic converters [4], [5]. 
However, to operating at high frequencies, it requires computational hardware must have equivalent 
processing capability. Previously, some studies have successfully implemented the field-oriented control 
algorithm using a DSP processor as shown in [6], [7]. This research focused on programming the pulse width 
modulation (PWM) pulse generator and current loop control. However, the higher the sampling frequency, 
the more computational resources are required, which leads to the time left to implement other functions, and 
control loops are limited [8]. FPGA technology [9] with the ability to handle multi-threading, which allows 
high frequency sampling and flexibility [10] is the solution for this problem. Furthermore, FPGA 
manufacturers also provide a system-on-chip (SoC) environment that making it increasingly used in power 
electronics and electrical drives systems [11]-[18]. In [19]-[21], the V/f control for induction motor using 
FPGA platform has been proposed. The advantage of this method is simple design, but V/f scalar control is 
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only suitable for low power and low control accuracy applications. The stator voltage regulator in [22] uses 
the sinusoidal pulse width modulation technique (SPWM) to generate gates pulse, which does not maximize 
the utilization of DC side voltage of three-phase voltage source inverter (VSI). This paper focuses on 
implementing the FOC algorithm using Kit FPGA Xilinx Zybo Zynq 7010 in VHDL language programmed 
by manual coding method for GaN-based induction motor drive. Figure 1 depicts the structure of the FOC 
algorithm that will be implemented in this work. 
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Figure 1. Block diagram of field-oriented control algorithm for induction motor drive 


2. IMPLEMENTING FOC ALGORITHM ON FPGA PLATFORM 

This section will present in detail the method of implementing the FOC algorithm on FPGA 
platform. All blocks in Figure 1 will be programmed and packaged as an integrated circuit (IC) in the 
development board. 


2.1. FOC algorithm and modeling the induction motor 

The FOC algorithm was developed decades ago by Siemens’s engineers [23]. The vector control 
algorithm is based on the idea that stator current can be decoupled into two components by using transformation 
from stationary abc frame (Isa, Isp, Isc) to rotating dq frame (Isa, Isq), the first one (Isa) produces flux and the other 
(Isq) produces torque. Thus, the induction motor now can be modeled as a separately excited DC motor and 
make DC motor control theories can be used for AC motor. Assume that the parameters of induction motor are 
stator resistance R,, stator lost inductance Los, mutual inductance Lm, rotor inductance R,, rotor lost inductance 
Los, number of pole pairs Zp, synchronous speed @s, rotor speed œ, rotor flux y, and T is sampling time. With 
LeLesthm, LeLatbm Os=@-0, T=L/T, T=LYTs ¢ =1—-12,/(LsL,), Wia = Wra/lm. 4 difference 
equations of induction motor in dq frame can be given as [24]: 
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In (1) and (2) are used to calculate the parameters of the current controller, where (3) and (4) are used to 
design the flux model of induction motor. 


2.2. Measurement of feedback current by XADC 


XADC is an ADC module integrated into some FPGA series of Xilinx, including a dual 12-bit, 
1 Mega sample per second (MSPS) ADC, and on-chip sensors [25]. Because induction motor is balance load so: 
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Isa+IsotIsc=0 (5) 


Therefore, only two input pairs are used to measure feedback current Isa and Is, then Ise will be 
calculated by using (5). In order to measure them, current signals must be transformed to voltage signals by: 


U(V) = 0.01 - I(A) + 0.5 (6) 


where I(A) is the instantaneous value of feedback current while U(V) is the value of voltage signal 
transformed. Since the ADC signal’s value is in the range of 0 to 4095, permissible current measurement 
range is from -50 to 50 A. 


2.3. Measuring the feedback rotor speed 

In this work, the rotor speed is measured by using the incremental encoder and will be decoded by 
the QEP resolver programmed in FPGA. Specifically, let consider the states of 2 channels A and B at ti, t2, ts, 
t4 in Figure 2. When the motor is in the forward or backward rotate condition, the states of t1, t2, t3, t4 are also 
listed respectively. When all 4 states are read, the encoder has rotated 1 pulse cycle, then combine 2 adjacent 
states to obtain the values as shown in Table 1. By reading these combined values, both the direction of 
rotation and rotor speed of the induction motor will be determined. 
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Figure 2. The logic state of 2 channels A and B 


Table 1. Combination of 2 adjacent states 


Combination state Forward rotation Backward rotation 
t-t 0100 1000 
t2- t; 0010 0001 
t3— t4 1011 0111 


t4-tı 1101 1110 


Now let use a variable to count encoder pulses named encodercount, another variable read the state 
of encoder named EncoderVal. When the EncoderVal’s value equal to any value in Table 1, the encodercount 
will be added or subtracted by 1 respectively. Otherwise, no changes are made to the variable encodercount. 
In this application, the clock frequency of this module is 1MHz and rotor speed will be sampled at 2000 Hz. 
Assume that when the counter is overflow, encodercount’s value is x, number of pole pairs is zp, the encoder 
has n PPR, the electrical rotor speed is calculated as: 


We = x - 20002, : 2/(4n) (rad/s) (7) 


2.4, Frame transformation 

As depicted in Figure 1, both the flux controller and the speed controller generate the reference for 
the current controllers in dq frame. They will be compared to the feedback current converted to the dq frame 
through park and clark transformation. The output of the current controller is Vsa and Vsq in dq frame will be 
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transformed to aß frame to go to the SVPWM module. The transformation between them is presented as the 


Table 2. The frame transformation modules are programmed in the FPGA platform by using (8) to (13). The 
trigonometric functions sin(@) and cos(0) will be determined by using the look-up-table method. 


Table 2. Fomula for frame transformation 


Clarke transform Park transform Invert park transform 
2 1 1 Isa = Isa ' Cos( 8) + Isg + Sin(0) (10) Usa = Usa» Cos( 0) — Ugg ' Sin( 0) (12) 
bsa = lsa'3 =la gt (© =h C i tu = Ugg Sin( 0) + Ugg * Cos(0) (13 
sq — 4sp* os( 0) —Isq + Sin(0) (11) sa = “sd in( )+ sq os( ) ( ) 
lsg = I.p/V3 = Isc/ V3 (9) 


2.5. Programming flux model 

Flux model the FOC algorithm is used to estimate the rotor flux, calculate the synchronous speed 
and the rotor angle. Firstly, rotor flux can be estimated directly from (3). The synchronous speed is calculated 
from (4): 


isq(k) 
w, (k) = w(k) + ae (14) 


Then the rotor angle is the result of integrating the synchronous speed œs by (15): 
O = T + ws-1) + O(K-1) (15) 
The difference (3), (14), and (15) can be implemented by a digital system. 


2.6. PI control algorithm 

In this project, the PI controller is applied for both inner and outer control loops, this is the most 
control algorithm and has been demonstrated to be effective for motor control. The PI controller is described 
in a differential as in (16): 


u(t) = fp(t) + fi) = Kp - e(t) + K; f e(t)dt (16) 


Where K, is the proportional gain, Kj is the integral time constant, e(t) is the steady-state error, and u(t) is the 
control signal of the controller. With the sample interval T, (16) also can become a difference equation by 
discretization as follows these steps: 


discretization 


With the proportional component: fp(t) = K, : e(t) ——————> fp(k) = Kp -e(k) (7) 


The discretization of integral component is implemented by following the Riemann Sums with the sample 
interval T. 


Ki 


filk) = 2+ Dkoe(n) = F e(k) +. Dkzbe(n) = F elk) + fk - 1) (18) 


T 


Finally, the difference equation of PI controller is calculated as: 
Ki 
u(k) = fi(k — 1) + Kpe(k) + > e(k) (19) 


From (19) can be applied for digital PI controller on any digital system such as FPGA, because PI 
controller is designed in the time domain, so the discretization can reduce the control quality of the 
controller. To avoid this phenomenon, the higher the sampling frequency, the more similar the characteristics 
of the controller designed in the digital domain to its characteristics in the time domain. FPGA platform is a 
powerful tool to achieve high control quality. 


2.7. SVPWM technique 

Space vector modulation (SVM) is a technique used for PWM, it is commonly applied for VSI. The 
biggest advantage of this technique is that it can utilize DC side voltage (~0.577 Uac) better than the SinPWM 
technique. In this research, it will be programmed to generate the gate pulse at frequency 100kHz for 2-level 
VSI using GaN devices. The SVPWM technique has been proposed in detail in [8]. Figure 3 depicts the 
SVPWM subsystem programmed in FPGA. Figure 4 depicts the flowcharts of the important blocks in this 
work to clarify the programming method, which presented in the previous sections. 
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Figure 3. Block diagram of SVPWM technique in FPGA 
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Figure 4. Flow charts of, (a) QEP resolver, (b) flux model, and (c) PI controller 


3. RESULTS AND DISCUSSION 
In this section, several of the results when launching the FOC algorithm implemented in the FPGA 


platform with Typhoon HIL 402 devices to verify the validity of this research and make some discussion. 
System clock of the design is 100MHz. The connecting diagram in detail of FPGA and Typhoon HIL is 
shown in Figure 5. The parameter is presented in Table 3 and Table 4. 


Table 3. Parameters of induction motor 
Rs Los Ln R: Lor Z J Pole pairs Ba 
2.52 Q 0.0062 H 0.1763 H 0.97 Q 0.0095 H 2 0.117 Kg.m/s? 2 1.2341 Wb 
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Table 4. Parameters of (a) VSI circuit and (b) controllers 


(a) (b) 
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Uae! 700 V Speed controller 5 178 +22 
Ci: 1000 uF Current controller 30 6690 +310 
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Figure 5. Connecting FPGA and typhoon HIL 


3.1. Device utilization and processing time 

Here, the device utilization of the FPGA board is counted in Table 5 giving an overview of the 
optimization of design. The processing time also calculated by setting a flag at the beginning and the end of 
the computation cycle, Figure 6 shown the processing time of FPGA in one sample. The processing time is 
0.22us, it significantly smaller than the sample time and allows the duty cycle value just calculated can be 
employed immediately after finishing the calculation instead of delay to the next cycle. 


Table 5. Device utilization in FPGA board 


Parameter Used Available Utilization 
Look up table (LUT) 9218 17600 53.8% 
LUTRAM 24 6000 0.4% 
Flip flop (FF) 2542 35200 7.2% 
BRAM 1 60 47.6% 
DSP 26 80 32.5% 
IO 17 100 17% 
BUFB 4 32 12.5% 
MMCM 1 2 50% 


3.2. Result when launching with HIL devices 

In this section, the magnetizing process of induction motor, the waveform of stator current, speed 
response, and motor mechanical torque when operating in 4 difirrence scenarios, magnetizing process, no- 
load condition, and constant load + 30Nm will be presented in Figures 7 to 11 to demonstrate the efficiency 
of FOC algorithm implemented on FPGA platform. 
Scenario 1: Magnetizing process 

Figure 7 shown the current response in magnetizing process, the magnetizing current is 7A as expected 
because it always equals to Isa and double times Is and Isc at the final value. The settling time is approximately 3 
ms and there is no overshoot and ringing. Figure 7 to 11 shown that the stator current when motor operating 
with different load has a sinusoidal waveform, rotor speed responses equivalent to the speed command as in 
theory in every scenario. The design has been confirmed that it meets the requirement of system. 
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Scenario 2: No-load condition 


Scenario 3: Operating with constant load 30Nm 
Scenario 4: Operating with constant load -30Nm 
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Figure 6. Procesing time in FPGA 
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Figure 7. Magnetizing process of induction motor 
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Figure 8. Stator current of 2”¢ scenario when speed command is (a) 190rad/s and (b) 20 rad/s 
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Figure 9. Stator current of 3" scenario when speed command is (a) 190rad/s and (b) 20 rad/s 
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Figure 10. Stator current of 4" scenario when speed command is (a) 190rad/s and (b) 20 rad/s 
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Figure 11. Speed response 


4. CONCLUSION 

The implementation of a vector control algorithm using FPGA for GaN-based induction motor drive 
was successfully developed by manual coding method in this contribution. From the achieved result, the 
processing time in one sample is 0.22us, which meets the requirement of the power electronics system with 
switching-frequency up to 100kHz. The magnetizing process is fast and has no overshoot and ringing, the 
induction motor can operate under different load conditions and different speed commands with only one 
setting for parameters of PI controllers. These results verify the efficiency of the algorithm programmed on 
the FPGA platform and become a worthy replacement for DSP or microcontroller platform. 
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